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A tunable double optomechanically induced transparency (OMIT) with a squeezed field is inves¬ 
tigated in a system consisting of an optomechanical cavity coupled to a charged nanomechanical 
resonator via Coulomb interaction. Such a double OMIT can be achieved by adjusting the strength 
of the Coulomb interaction, and observed even with a single-photon squeezed field at finite temper¬ 
ature. Since it is robust against the cavity decay, but very sensitive to some parameters, such as the 
environmental temperature, the model under our consideration can be applied as a quantum ther¬ 
mometer for precision measurement of the environmental temperature within the reach of current 
techniques. 

PACS numbers: 42.50.Wk, 46.80.-|-j, 41.20.Cv 


I. INTRODUCTION 

Electromagnetically induced transparency (EIT) is a 
kind of effect with a narrow transparency window within 
an absorption line of atoms ^ , due to quantum interfer¬ 
ence between two quantum pathways in A-type atoms. 
This effect plays a key role in modern quantum optics 
experiments and applications, such as enhanced nonlin¬ 
ear susceptibility [2], optical switch [3], slow and fast 
lights [4], quantum memory quantum interference 

[5] and vibrational cooling [^. Recently, the study of 
the EIT has been extended to multi-channels, e.g., the 
double EIT [lOl [TT] , and focused on simulation of other 
physical phenomena, including Anderson localization [12] 
and quasi-charged particles m- 

The EIT analog occurring in an optomechanical system 
is called optomechanically induced transparency (OMIT) 
which is caused by quantum interference between 
two quantum channels in a A-type hybrid level config¬ 
uration composed of photon states of the cavity and 
phonon states of the optomechanical resonator da US]. 
The OMIT has been explored both theoretically and 
experimentally Similar to the EIT, the study 

of the originally defined OMIT [2T] was extended to the 
double OMIT by two coupled optomechanical resonators 
[33] or by an optomechanical resonator coupled to other 
systems |331|31]. Besides, the double OMIT was explored 
from the fixed double OMIT [331 ES] [35] to the tunable 
one involving a controllable coupling [33] . 

The present paper intends to investigate the unique 
behavior of a double OMIT with a squeezed field. Due 
to involvement of the squeezed field, the OMIT is ro¬ 
bust against quantum noise and thus possible to be a 
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candidate of quantum memory |26| . But if the model is 
extended to be a double OMIT in a tunable manner, the 
physics turns to be largely different. The key point is 
that the double OMIT is robust to the cavity decay and 
quantum noise of the environment can be correlated to 
the temperature-dependent noise. As such, we may carry 
out precision measurement of the environmental temper¬ 
ature, assisted by the squeezed field and the homodyne 
spectroscopy. Using other unique characteristics, preci¬ 
sion measurements of other parameters of the system are 
also available. This implies that the model under our 
consideration is by no means a simple extension of the 
previously considered OMIT [26], but with much differ¬ 
ent characteristics and applications. 

The temperature dependence is from the quantum field 
involved in our double OMIT, by which we are able to 
know the environmental temperature through detecting 
the noise spectra of the optomechanics m- This is very 
different from the OMIT with classical lights [HIMIIIH], 
whose noise spectra have no relevance to the environmen¬ 
tal temperature even under the cryogenic condition [14] . 
In this context, our scheme is also very different from the 
previous OMIT measurements resulted from the prop¬ 
erties of the OMIT spectra im [31]. As a result, our 
scheme provides a new paradigm for precision measure¬ 
ment based on the noise imiM]. On the other hand, 
compared with the conventional OMIT |3^, our double 
OMIT characterizes as a linear variation of the peak value 
with respect to the environmental temperature, which is 
robust to the cavity decay and does not vary with the 
Coulomb coupling between the two nano-mechanical res¬ 
onators (NAMRs). This feature exists no matter whether 
the two NAMRs are identical or not, which is useful for 
practical applications as discussed later. 

The paper is organized as follows. In Sec. II, we 
present the solution to the model of our designed dou¬ 
ble OMIT and focus on the spectra via a homodyne de¬ 
tection. In Sec. Ill, some calculations are made nu- 
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merically with experimentally available values, justify¬ 
ing some unique features, such as robustness against the 
cavity decay and invariance with the Coulomb coupling 
strength. The feasibility of precision measurement of the 
environmental temperature is discussed in Sec. IV. Other 
extended discussion is made in Sec.V and a brief conclu¬ 
sion is given in the last section. 


II. THE MODEL AND SOLUTION 




FIG. 1: (Color online) Schematic diagram of the double 
OMIT system and the measurement. A high-quality Fabry- 
Perdt cavity consists of two fixed mirrors and a charged 
NAMRi, which is charged by the bias gate voltage Vi and 
subject to the Coulomb force due to the charged NAMR 2 
outside the cavity and with the bias gate voltage V 2 . The op¬ 
tomechanical cavity of length L is driven by two light fields, 
one of which is the pump field ei with frequency uii and the 
other of which is the probe field Ci„ with frequency ujp. qi and 
52 represent, respectively, the small displacements of the two 
NAMRs from their equilibrium positions, with ro the equilib¬ 
rium distance between them. The output field Cout from the 
cavity turns to be Cout , which is mixed with a strong local field 
cio centered around the probe frequency Up at a 50:50 beam 
splitter (BS). Finally the homodyne spectra are obtained by 
the spectrum analyzer (SA) assisted by two photon detectors 
(PD). 


As sketched in Fig. 1, there are two charged NAMRs 
with one (NAMRi) inside a Fabry-Perdt (FP) cavity cou¬ 
pling to the cavity mode by the radiation pressure and 
interacting with the other (NAMR 2 ) outside the cavity. 
The FP cavity contains two mirrors distant by L with 
the left-hand side mirror partially transmitting and the 
right-hand side one 100% reflecting. There is a driving 
on the cavity mode from the left-hand side mirror by a 
strong coupling field with frequency wj. The system can 


be described as 

Hwhole — ?iCUcC^C-l- ^ ^ 2jn 

i=i 

— hgc^cqi + Hi + ihei{c^ — H.C.), 

where the first term is for the single-mode cavity field 
with frequency U3c and annihilation (creation) operator 
c(c^). The second and third terms describe the vibration 
of the charged NAMRs with frequency 0 Ji(u} 2 ) and effec¬ 
tive mass toi(to 2 ). Pi{p 2 ) and qi(q 2 ) are the momentum 
and the position operators of NAMRi (NAMR 2 ), respec¬ 
tively. The fourth term presents the radiation pressure 
coupling the cavity field to the NAMRi with a coupling 
strength g = ujdL. 

Coulomb coupling between the two charged NAMRs is 
given by Hi = , where ro is the distance 

between the equilibrium positions, and NAMRi and 
NAMR 2 take the charges CiVi and —C 2 V 2 with Ci and 
C 2 being the capacitance of the gates, respectively. Un¬ 
der the assumption that the deformations of the NAMRs 
are much less than their distance ( 51,(72 ^ ?'o)i the 
Hamiltonian Hi can be expanded to the second order as 
Hi = - ^ + (“)"]; Since the linear 

term may be absorbed into the definition of the equi¬ 
librium positions, and the quadratic term includes the 
renormalized oscillation frequencies for both the NAMRs, 
we have a reduced form Hi = hXqiq 2 for A = 

EIllSl]. 

The last term represents the cavity field driven by 
an input field with frequency wi, where the pump field 
strength e; = depends on the power p of the 

coupling field and the cavity decay rate k. 

In a frame rotating with the pump field frequency wp 
the Hamiltonian in Eq. ([^ is rewritten as 

Htotai = HAcc'^c-b ^(^^ + 

1=1 ^"^3 

— figc^ cqi + fi\qiq 2 + ih£i{c^ — c), (2) 


with Ac = LOc — oJi- Considering the decay rates 71 and 72 
for the NAMRi and NAMR 2 , respectively, we obtain the 
corresponding frequency-domain correlation functions for 
the thermal noise ^1 and ^2 at a temperature T, 


(Cr(w)^r(II)) = 2TTnjrmrU} 


1 + coth 


r nuj \ 
Wbt) 


^(w-l-H), 


where r = 1, 2 and fcs is the Boltzmann constant. 

We assume that the cavity mode c couples to the input 
quantum field Cm, which is a narrow-band squeezed field 
with the center around the frequency ojp = ujc + Wi and 
with a finite bandwidth F. The nonvanishing correlation 
functions for this input squeezed field are given by 

MF^ 

(c™(a;)ci„(H)) = 27r - ^6{oj + il. - 2 uji), 

I -I- (w - Wi)"^ 
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(ci„(a;)4„(-0)) = 2 tt 


NT^ 


r2 + (w-wi)2 


+ 1 


5 (io + ri), 


( 3 ) 


where N is the photon number in the squeezed vacuum, 
and M = N{N + 1 ) is an anti-normally ordered term 
including a broadband contribution from the vacuum 
noise. 

Considering the input squeezed field, Eq. ([^ un¬ 
der dissipation and fluctuation is governed by quantum 
Langevin equations, yielding 


qi 

c 

Pi 

P2 


Pi . P2 

mi m 2 

-[k + i(Ac - gqi)]c + ei + V^Cin, 

-niiujfqi - hXq2 + Ugc^c - yipi -I- ^i, 

-m 2 Ujlq 2 - hXqi - 72P2 + ^2, ( 4 ) 


whose steady solutions are given by. 


Pis 


P2s = 0 , 


qis = 


%|CsP 


miiJl — 
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with A — Ac — gqis being the effective detuning between 
the cavity and the driven fields. 

Next,we consider the linear operators as steady mean 
values plus additional fluctuation operators, 

qT = qrs+Sqr, PT=PTS+Spr, C = Cs + 6 c, 

where Sqr, Spr, Sc are small fluctuations around the cor¬ 
responding steady values. Since the steady values have 
no contribute on the output fields, we only focus on the 
fluctuation operators, which work as the probe fields and 
influence the output fields. 

Defining the fluctuation operators X = 
{ 6 pi{u}), Sqi{u}),Sp2{uj), Sq2{w), 5 c(a;), ( 5 c'l'(a;))'^, we 

reach the linearized quantum Langevin equations 
AX = B for the fluctuation operators from Eq. ( 4 ), 
where 
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and B = ( 0 , .^1(0;), 0 ,6(w), v^cj„(-w))^. 

Thus the fluctuation < 5 c(a;) of the cavity field can be 
solved by the linearized equations above. 

Based on the input-output relation Coutipj) = 
x/^c{ijj) — Cin{oj), we define the output field as Cout{^) = 
Coutipj) + Cin{oj) = x/^c{u}) in order to study the physi¬ 
cal properties of the light field leaking out of the cavity, 
as in [niEH!. Straightforward deduction yields 

6cout{uj) = E{uj)ciniuj) + F{uj)cy{-uj) 

+ l^l(w)^l(w) -b V'2(w)C2(w), (7) 


where 


E{ui) 


2 k[ - 77 --r 

''K + i{A — uj) 

^ifig‘^\csy [A + ui + iK)m 2 B 

(A — UJ — iK)d(uj) 


( 8 ) 


E(uj) 


ViM 

V2(uj) 


—2iKhg^Cgm2Bi 

d{uj) ’ 

x/^iCsg[—K -b i{A -b a;)]m2Bi 
d{uj) 

—^ 2 KhgCsX(A -b w -b in) 
d(uj) 


( 9 ) 


with d{u:) = —‘h?X^A + m 2 Bi{ 2 \cs\^g^fiA~\-miAB 2 ), A = 
A^-b(K —Bi = a;^-b*W72—and B2 = w^-b*a;7i — 

Based on the above mentioned correlation functions 
of Cin{uj) and Cr(w) as well as the standard homodyne 
detection [ 53 ] as plotted in Fig. 1 , we may understand 
characteristics of the system from the homodyne spec¬ 
trum X{uj), which can be analytically expressed as be¬ 
low provided that the fast oscillating terms at frequencies 
±2a;i(2) are omitted. 


X{uj) 
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E*{-oj + oj,)E*ioj + a;i)-—^ + \E{-uj + 


\E{uj + wi)p + \F{-uj + wi)|^ 

IT/ /.™ / .\ , 


11 - 1 - coth 

^(u! + UJi) 

‘Ih^T 

11 - 1 - coth 

Tl(bJ -\- Wi) 

[ 2kBT \ 


\Vi{—oj + oJi)\^Ti"fimi{oj — wi) X |l + coth 
\V2(—U! + wi)|^?i72m2(w — wi) X < 1 + coth 


h(uj — uji) 
2kBT 
h(ui — uji) 
2kBT 


( 10 ) 


where the first four terms are from the input squeezed 
field and the next two terms are relevant to the sponta¬ 
neous emission of the input vacuum noise. The rest terms 
are caused by the thermal noise of the NAMRs, which is 
temperature dependent. So the squeezed field employed 
does not work better for the temperature-dependent ef¬ 
fects in comparison to other quantum fields, but enhanc¬ 
ing the measurement precision in the homodyne spec¬ 
trum. In addition, Eq. (10) is more general with respect 
to the counterpart in Ref. [55], since it can be reduced 
to Eq. (13) in [^ if A = 0, i.e., in the absence of the 
Coulomb coupling. 


III. HOMODYNE SPECTRA OF THE DOUBLE 
OMIT 



We specify below some unique characteristics of the 
double OMIT by the numerically calculated homodyne 
spectra. For simplicity, we first suppose the two NAMRs 
to be identical in our treatment. The non-identical case, 
which is more general but not fundamentally different, 
will be justified later. 

Our numerical calculation is carried out based on re¬ 
alistic parameter values [34]. We consider an optome¬ 
chanical cavity with length L = 25 mm and decay rate 
K ~ 27r X 215 kHz, driven by the pump field of wave¬ 
length A ; = 2ttc/uji = 1064 nm. For the two identical 
NAMRs, we assume the effective mass mi = m 2 = 145 
ng, the eigen-frequencies Wm = wi = 022 = 27r x 947 kHz, 
the decay rates 7 ™ = 71 = 72 = 27r x 141 Hz, and the 
quality factors Qi = Q 2 = = 6700. In addition, 

the linewidth of the squeezed vacuum is supposed to be 

r = 2 k. 

In most of the calculations below, we employ the zero 
temperature T = 0 and the photon number N = 5 
in the squeezed vacuum, and assume the coupling field 
power p = 2 mW and the coupling strength unit Aq = 
4 X 10^® Hz/m^. The homodyne spectrum X{w) plotted 
in Fig. ia) presents the change from a single trans¬ 
parency window to two transparency windows with in¬ 
creasing Coulomb coupling, which reflects a fact that the 


FIG. 2: (Color online) The homodyne spectra X{uj) of the 
output field as functions of the normalized frequency ujjuirn, 
where (a) for different Coulomb coupling A with the coupling 
strength unit Ao = 4 x 10^® Hz/m^; (b) for different temper¬ 
ature T with N=5; (c) for different photon number N with 
T = IQ mK; (d) for different cavity decay k with the decay 
unit Ko = Stt X 215 kHz. 


Coulomb coupling breaks down the original interference 
in the OMIT and splits the bosonic mode of the system 
into two. Since the energy difference between the two 
split bosonic modes depends on the Coulomb coupling, 
the splitting of the transparency windows is relevant to 
the Coulomb coupling (53]. However, the middle peak is 
fixed no matter how much the Coulomb coupling varies. 

For a finite temperature, the homodyne spectrum 
X{uj) still works for the double OMIT, as presented in 
Fig. 2(b) where the visible middle peak and two nadirs 
exist even at T = 100 mK. However, the trend reflected in 
Fig. 2(b) indicates that the double-OMIT will definitely 
disappear with further increase of the temperature. Be¬ 
sides, a more careful calculation shows that the photon 
number plays an important role in the X(u}) variation. 
From Fig. 2(c), we find that the middle peak and two 
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nadirs are visible at T = 10 mK, even for the squeezed 
state at the single-photon level (the red dashed curve). In 
particular, compared with Fig. 2(b), we find that the two 
nadirs of the double OMIT are fixed with the variation of 
the photon number, but changing with the temperature. 

Moreover, the cavity decay modifies the profiles of the 
transparency windows, as plotted in Fig. id). However, 
although the profiles become narrower and sharper with 
smaller cavity decay rate k, the peak and the nadirs of the 
two transparency windows remain unchanged, implying 
robustness against the cavity decay at these points. 



FIG. 3: (Color online) (a) The homodyne spectra X{uj) as 
functions of the frequency uj/uJm and the temperature T for 
A = Aq. (b) The rescaled middle peak value Xp (in units 
of Xo = 87.59) as a function of the temperature T for A = 
Ao, where the temperature measurement is available within 
the range 0 < T <0.13 K. The upper limit of the measured 
temperature is restricted by the measurement resolution 1% 
of Xp at T =0.13 K. Other parameters take the same values 
as in Fig. 2(a). 


IV. MEASUREMENT OF THE 
ENVIRONMENTAL TEMPERATURE 

Using above indicated properties of the homodyne 
spectrum X{ui), we may carry out precision measure¬ 
ment of the environmental temperature using the double- 
OMIT with the squeeze field. 

Fig. [^a) presents the step-by-step change of X{w) 
with the environmental temperature, where the mid¬ 
dle peak increases linearly with the temperature. Since 
we have renormalized the middle peak value by Xp = 
A(0)/Xo with Xq = 87.59 = A(0) at zero temperature, 
the temperature change can be exactly known from Fig. 
3(b) by precisely measuring the variation of Xp. But as 
shown in Fig. 3(a), the height difference between the 
middle peak and the two nadirs shrinks with the tem¬ 


perature increasing. This implies an upper limit of the 
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FIG. 4: (Color online) The quantum signal visibility {QSV) 
as functions of the photon number N and temperature T 
for A = Aq. The dashed line in green represents a border¬ 
line for available measurement, below which QSV is less than 
the measurement resolution 1%. Other parameters are of the 
same values as in Fig. 2(a). 


measured temperature, e.g., T = 0.13 K where the dou¬ 
ble OMIT reaches the resolution limit of the observation. 

Specifically, for the linear variation of Xp with respect 
to T in Fig. ib), the sensitivity can be evaluated by the 
slope k = dT/dX{0) = 4.4 x 10“^K. As a result, for a 
measurement resolution dA(0) = 1% of the peak value, 
the detectable temperature change can be minimized to 
4.4 X 10“® K, which is lower by one order of magnitude 
compared to a quantum thermometer designed based on 
a noise measurement using the electron charge [35]. The 
linear variation of the middle peak with respect to the 
temperature can be understood from Eq. (10) in which 
the last four terms are reduced to be linearly changing 
with T if T —>-0. In fact, the measurement precision in 
our case can be further enhanced if we elaborately change 
the mass ratio of the two NAMRs, as discussed later. 
We have to emphasize that this measurement based on 
the peak values is insensitive to the change of Coulomb 
coupling, as indicated in Fig. 2(a). 

To carry out a precision measurement in our scheme, 
we have to have a big enough contrast of X{uj) for our 
observation. To this end, we consider below the influence 
from the photon number and the environmental temper¬ 
ature. When the two NAMRs are identical, if many pho¬ 
tons are involved in the squeezed state, we have N « M, 
and Eq. (10) at low temperature (T —>• 0) expands to the 
first order at the frequency of the peak point {lOpeak = 0) 
as 
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X{0Jpeak)=X{0) = N[E{0Jm)+E*iu^)]^ + \EM\^ + \F{0Jm)\^ 

+ 2\Vi{uJm)\‘^fl^imiUJmC0th[ ] 

+ 2|l/2(wm)P?i72TO2a;mC0th[^^^] 

= N[E(u}rn) + E* {uJm)]'^ + \E(uJm)\'^ + \E{LiJm)\‘^ + k{0)T, (11) 


where the slope is fc(0) = 4:kBjm'm[\Vi{u}m)\‘^ + 

|h 2 (wm)P], N is relevant to the quantum signal terms 
which compete with the last two thermal noise terms in¬ 
volving T. With the increase of T, the values of the 
thermal noise terms will exceed those of the quantum 
signal terms and thus the quantum signal is completely 
buried by the thermal noise, i.e., disappearance of the 
double OMIT. To clarify this point, we define quantum 
signal visibility QSV = (peak — nadir)/peak as a con¬ 
trast of our observation, where QSV = 1 implies an 
overwhelming quantum signal and disappearance of the 
double OMIT corresponds to QSV = 0. Involvement of 
more photons helps increasing the contrast, as indicated 
in Fig. 2(c) and understood from Eq. (11). We plot in 
Fig. 4 a borderline for available precision measurement, 
below which QSV is smaller than 1%, the measurement 
resolution we assumed above. As such, if T = 0, the 
double OMIT always exists no matter how many pho¬ 
tons are involved. But with T increasing, more photons 
are required to be involved if our proposed measurement 
of the environmental temperature works. 

V. DISCUSSION 

The precision measurement of the environmental tem¬ 
perature described above is based on identical NAMRs. 
A more general and realistic case is with non-identical 
NAMRs, which are different in frequency or mass. In 
such cases, our model has different characteristics and 
thus different applications. 

As an example, we hrst consider in Fig. 5 the situ¬ 
ation with different frequencies of the two NAMRs. In 
this case, the profile of the double OMIT keeps changing 
with the frequency difference Auj = |wi — UI 2 \^ where the 
double windows are first split into triple windows and 
then become a standard OMIT with a single window, 
like the absence of Coulomb coupling (See Fig. 5(a)). 
Fig. 5(b) provides another view angle to observe the 
role played by the Coulomb coupling. Different from 
the situation with identical NAMRs, the Coulomb cou¬ 
pling yields triple windows directly from the single win¬ 
dow once the Coulomb coupling turns to be non-zero. 
This feature is actually resulted by the interference of 
two double OMITs with two asymmetric windows due 
to the frequency difference. As indicated by the homo¬ 
dyne spectra in Figs. 5(a) and 5(b), although the profiles 
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FIG. 5: (Color online) (a) The homodyne spectra X{ljj) as 
functions of oj/lji and Ata/cai for A = Ao- (b) The homodyne 
spectra A(a;) as functions of oi/oJi and A/Ao with Auu — O.lcJi. 
(c) The rescaled peak value Xp (in units of 84.03) as a function 
of the temperature T for A = Aq, where the available mea¬ 
surement of the temperature is 0 < T <0.14 K (Aca = O.lwi) 
or 0 < r <0.18 K {Auj = 0.2tai). The upper limit of the 
measured temperature is restricted by the measurement res¬ 
olution 1% of Xp. Other parameters are of the same values 
as in Fig. 2(a). 

of the spectra change in the variation of Aw and A, the 
two symmetric peaks in the case of triple windows are al¬ 
ways fixed, where the two peak values can be evaluated 

by ^^\ui=uj± = 0 and < q, with w± the 

frequencies relevant to the two symmetric peaks. 

Based on this feature, we consider below a measure¬ 
ment of the environmental temperature using one of the 
peak values (See Fig. 5(c)). The peak value varies lin¬ 
early with respect to T, in the same fashion as in Fig. 
3(b) for the middle peak values in the case of identical 
NAMRs. We may find a similar formula to Eq. (11) for 
the slope of the peak values X{ujpeak) varying with the 
environmental temperature T as 

2 

k{uJpeak) — ‘2kB ^ ^ Tj '^j [ \ Vj {^peak T ^l)| 

i=i 

4”|V/( Wpeafc 4” Wi)| ] = k{ ^peak)- (1^) 

Straightforward calculations of the slopes in this case 
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FIG. 6: (Color online) (a) The homodyne spectra X{uj) as 
functions of ujjuJm for A = Ao and T=0. (b) The homodyne 
spectra X{uj) as functions of ujjuJm for A = Ao and T = 20 
mK. (c) The rescaled middle peak value Xp (in units of 87.59) 
as a function of the temperature T for A = Ao, where the 
temperature measurement is available within the range 0 < 
T <0.16 K (m 2 = 0.5mi) or 0 < T <0.12 K (m 2 = 2mi). 
The upper limit of the measured temperature is restricted by 
the measurement resolution 1% of Xp. Other parameters take 
the same values as in Fig. 2(a). 

present less precise measurements of the environmental 
temperature compared to the case of identical NAMRs 
since we have the sensitivity with 0.44 (Aw = O.lwi) 
or 0.33 (Aw = 0.2wi) of the counterpart in the case 
of Aw = 0. Therefore, for a more precise measure¬ 
ment of the environmental temperature with two differ¬ 
ent NAMRs, the frequency difference is required to be as 
tiny as possible. 

If the two different NAMRs are with the same fre¬ 
quency but with different mass, we have only double 
OMITs, rather than triple OMITs. In this case, we found 
that the middle peak remains the same value for differ¬ 
ent ratios of m 2 /mi provided that the temperature is 
zero, but varies with different slopes for different ratios 
of m 2 /mi if T 7^0 (See Fig. 6 ). In particular, for a bigger 
mass of NAMR 2 than NAMRi, e.g., m 2 = 2mi in Fig. 
6 (c), the measurement sensitivity of the environmental 
temperature is higher than the counterpart in the iden¬ 


tical case. In addition, the measurement range changes 
with different ratios of TO 2 /TO 1 as indicated in Fig. 6 (c). 
Nevertheless, within the range T G [0,0.12] K, we have 
the possibility to obtain the measurement precision of 
temperature better than 4.4 x 10“® K. 

VI. CONCLUSION 

In summary, we have justified the possibility of pre¬ 
cisely detecting the environmental temperature by the 
unique quantum characteristics of double-OMIT. To our 
knowledge, this is the first scheme for such a precision 
measurement in the optomechanical system via the noise 
spectra. We have shown by numerical simulation that 
we are able to have a better precision of measuring the 
environmental temperature than a previously reported 
quantum thermometer |35j . 

For simplicity, however, we have remained the radia¬ 
tion pressure as a constant throughout the paper. For a 
thorough investigation of the temperature measurement, 
it is necessary to explore the change of the radiation 
pressure. Straightforward calculations indicate that en¬ 
hancement of the radiation pressure due to increase of 
the photon number will definitely lead to a better pre¬ 
cision measurement of the environmental temperature. 
Nevertheless, the nonlinear effect in the optomechanics 
due to more photons involved would bring in unexpected 
complexity, which needs further scrutiny. 

Moreover, due to the tunable fashion and robustness to 
cavity decay, the model under consideration can also be 
applied to other applications, such as precisely measur¬ 
ing Coulomb coupling strength and the frequency (mass) 
difference between the two NAMRs. Further exploration 
would be more interesting and is underway. 
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